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When we are fortunate enough to view an exoplanetary system nearly edge-on, the star and planet 
periodically eclipse each other. Observations of eclipses — transits and occultations — provide a bonanza 
of information that cannot be obtained from radial-velocity data alone, such as the relative dimensions 
of the planet and its host star, as well as the orientation of the planet's orbit relative to the sky plane and 
relative to the stellar rotation axis. The wavelength-dependence of the eclipse signal gives clues about the 
the temperature and composition of the planetary atmosphere. Anomalies in the timing or other properties 
of the eclipses may betray the presence of additional planets or moons. Searching for eclipses is also a 
productive means of discovering new planets. This chapter reviews the basic geometry and physics of 
eclipses, and summarizes the knowledge that has been gained through eclipse observations, as well as the 
information that might be gained in the future. 
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1. INTRODUCTION 

Frotn immemorial antiquity, men have dreamed of a 
royal road to success — leading directly and easily to some 
goal that could be reached otherwise only by long ap- 
proaches and with weary toil. Times beyond number, this 
dream has proved to be a delusion.... Nevertheless, there 
are ways of approach to unknown territory which lead sur- 
prisingly far, and repay their followers richly. There is 
probably no better example of this than eclipses of heav- 
enly bodies. — Henry Norris Russell (1948) 

Vast expanses of scientific territory have been traversed 
by exploiting the occasions when one astronomical body 
shadows another The timing of the eclipses of Jupiter's 
moons gave the first accurate measure of the speed of light. 
Observing the passage of Venus across the disk of the Sun 
provided a highly refined estimate of the astronomical unit. 
Studying solar eclipses led to the discovery of helium, the 
recognition that Earth's rotation is slowing down due to 
tides, and the confirmation of Einstein's prediction for the 
gravitational deflection of light. The analysis of eclipsing 
binary stars — the subject Russell had in mind — enabled a 
precise understanding of stellar structure and evolution. 

Continuing in this tradition, eclipses are the "royal road" 
of exoplanetary science. We can learn intimate details about 
exoplanets and their parent stars through observations of 
their combined light, without the weary toil of spatially re- 
solving the planet and the star (see Figure [T]i. This chapter 
shows how eclipse observations are used to gain knowledge 
of the planet's orbit, mass, radius, temperature, and atmo- 
spheric constituents, along with other details that are other- 
wise hidden. This knowledge, in turn, gives clues about the 
processes of planet formation and evolution and provides a 
larger context for understanding the properties of the solar 
system. 

An eclipse is the obscuration of one celestial body by an- 
other When the bodies have very unequal sizes, the passage 
of the smaller body in front of the larger body is a transit 



and the passage of the smaller body behind the larger body 
is an occultation. Formally, transits are cases when the full 
disk of the smaller body passes completely within that of 
the larger body, and occultations refer to the complete con- 
cealment of the smaller body. We will allow those terms to 
include the grazing cases in which the bodies' silhouettes 
do not overlap completely. Please be aware that the exo- 
planet literature often refers to occultations as secondary 
eclipses (a more general term that does not connote an ex- 
treme size ratio), or by the neologisms "secondary transit" 
and "anti-transit." 

This chapter is organized as follows. Section 2 de- 
scribes the geometry of eclipses and provides the founda- 
tional equations, building on the discussion of Keplerian or- 
bits in the chapter by Murray and Correia. Readers seeking 
a more elementary treatment involving only circular orbits 
may prefer to start by reading Sackett (1999). Section 3 dis- 
cusses many scientific applications of eclipse data, includ- 
ing the determination of the mass and radius of the planet. 
Section 4 is a primer on observing the apparent decline in 
stellar brightness during eclipses (the photometric signal). 
Section 5 reviews some recent scientific accomplishments, 
and Section 6 offers some thoughts on future prospects. 

2. ECLIPSE BASICS 

2.1 Geometry of eclipses 

Consider a planet of radius Rp and mass Mp orbiting a 
star of radius i?* and mass Af*. The ratio Rp/R^, occurs 
frequently enough to deserve its own symbol, for which we 
will use k, in deference to the literature on eclipsing binary 
stars. As in the chapter by Murray and Correia, we choose 
a coordinate system centered on the star, with the sky in 
the X-Y plane and the +Z axis pointing at the observer 
(see Figure |2]i. Since the orientation of the line of nodes 
relative to celestial north (or any other externally defined 
axis) is usually unknown and of limited interest, we might 
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star - planet shadow 

Fig. 1. — Illustration of transits and occultations. Only the combined flux of the star and planet is observed. During a transit, the flux 
drops because the planet blocks a fraction of the starlight. Then the flux rises as the planet's dayside comes into view. The flux drops 
again when the planet is occulted by the star. 



as well align the X axis with the line of nodes; we place the 
descending node of the planet's orbit along the +X axis, 
giving n = 180°. 

The distance between the star and planet is given by 
equation (20) of the chapter by Murray and Correia: 



a(l ~e^) 
1 + e cos / ' 



(1) 



where a is the semimajor axis of the relative orbit and / 
is the true anomaly, an implicit function of time depending 
on the orbital eccentricity e and period P (see Section 3 of 
the chapter by Murray and Correia). This can be resolved 
into Cartesian coordinates using equations (53-55) of the 
chapter by Murray and Correia, with ft = 180°: 



X = — rcos(a; + /), 

Y — — rsin(LLi + /) cosi, 

Z = r sin(aj + /) sini. 



(2) 
(3) 
(4) 



If eclipses occur, they do so when rgj^y = \/ X^ + is 
a local minimum. Using equations (|2][3]l. 



''sky 



1 + e cos / V 



(5) 



Minimizing this expression leads to lengthy algebra (Kip- 
ping 2008). However, an excellent approximation that we 
will use throughout this chapter is that eclipses are centered 



around conjunctions, which are defined by the condition 
X ^ and may be inferior (planet in front) or superior 
(star in front). This gives 



/tra ^" 2 /occ — 



2 



(6) 



where here and elsewhere in this chapter, "tra" refers to 
transits and "occ" to occultations. This approximation is 
valid for all cases except extremely eccentric and close-in 
orbits with grazing eclipses. 

The impact parameter b is the sky-projected distance at 
conjunction, in units of the stellar radius: 

6tra = ( — -. I , (7) 



a cos I 
acosi 



1 + e sin a; 
1 — e sin LJ 



(8) 



For the common case i?* <C a, the planet's path across 
(or behind) the stellar disk is approximately a straight line 
between the points X = ±R^,^/T~^ at F = 

2.2 Probability of eclipses 

Eclipses are seen only by privileged observers who view 
a planet's orbit nearly edge-on. As the planet orbits its star, 
its shadow describes a cone that sweeps out a band on the 
celestial sphere, as illustrated in Figure [3] A distant ob- 
server within the shadow band will see transits. The open- 
ing angle of the cone, Q, satisfies the condition sinG = 
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Fig. 2. — Illustration of a transit, showing the coordinate system discussed in Section 2.1, the four contact points and the quantities T 
and T defined in Section 2.3, and the idealized light curve discussed in Section 2.4. 



(i?* + Rp)/r where r is the instantaneous star-planet dis- 
tance. This cone is called the penumbra. There is also an in- 
terior cone, the antumbra, defined by sin Q — [R^, — Rp)/r, 
inside of which the transits are full (non-grazing). 

A common situation is that e and lu are known and i is 
unknown, as when a planet is discovered via the Doppler 
method (see chapter by Lovis and Fischer) but no informa- 
tion is available about eclipses. With reference to Figure |3] 
the observer's celestial longitude is specified by uj, but the 
latitude is unknown. The transit probability is calculated as 
the shadowed fraction of the line of longitude, or more sim- 
ply from the requirement \b\ < 1 + k, using equations (|7][8]l 
and the knowledge that cos i is uniformly distributed for a 
randomly-placed observer. Similar logic applies to occulta- 
tions, leading to the results 

fR^±Rp\ /l + esina;\ 

where the "+" sign allows grazing eclipses and the "— " 
sign excludes them. It is worth committing to memory the 
results for the limiting case Rp <C -R* and e = 0: 

p.,,.p„„.5l.O,005 (|l)(j|jj)-\ ,11) 



For a circular orbit, transits and occultations always go to- 
gether, but for an eccentric orbit it is possible to see transits 
without occultations or vice versa. 

In other situations, one may want to marginalize over 
all possible values of oj, as when forecasting the expected 
number of transiting planets to be found in a survey (see 
Section 4.1) or other statistical calculations. Here, one can 
calculate the solid angle of the entire shadow band and di- 
vide by Att, or average equations (|9lfT0ll over cj, giving 

Suppose you want to find a transiting planet at a particu- 
lar orbital distance around a star of a given radius. If a frac- 
tion rj of stars have such planets, you must search at least 
N {ri ptra)~^ Stars before expecting to find a transiting 
planet. A sample of >200 77"^ Sun-like stars is needed to 
find a transiting planet at 1 AU. Close-in giant planets have 
an orbital distance of approximately 0.05 AU and rj « 0.01, 
giving >10^ stars. In practice, many other factors affect 
the survey requirements, such as measurement precision, 
time sampling, and the need for spectroscopic follow-up 
observations (see Section 4.1). 

2.3 Duration of eclipses 

In a non-grazing eclipse, the stellar and planetary disks 
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Fig. 3. — Calculation of the transit probability. Left. — Transits are visible by observers within the penumbra of the planet, a cone with 
opening angle O with sin O = (7?* + Rp)/r, where r is the instantaneous star-planet distance. Right. — Close-up showing the penumbra 
(thick lines) as well as the antumbra (thin lines) within which the transits are full, as opposed to grazing. 



are tangent at four contact times tj-tiv, illustrated in Fig- 
ure |2] (In a grazing eclipse, second and third contact do 
not occur) The total duration is Ttot = ^iv — ^i, the 
full duration is Tf^w = tin — in, the ingress duration is 
■'"ing = tu — ti, and the egress duration is Tegr = tiy — tm. 

Given a set of orbital parameters, the various eclipse du- 
rations can be calculated by setting equation (|5]l equal to 
i?* ± Rp to find the true anomaly at the times of contact, 
and then integrating equation (44) of the chapter by Murray 
and Correia, e.g.. 



p /■/ill 




27rVl - e2 Jf,, 


a 



df. 



For a circular orbit, some useful results are 



P 



R, v/(l + fc)2-62 



(13) 



(14) 



Tfuii — hu — tu — — sin 

TT 



R^^iT^W^ 



(15) 

For eccentric orbits, good approximations are obtained by 
multiplying equations ( 11411151 ) by 



Xif,) [e = 0] _ 



lie sin Ll) ' 



(16) 



a dknensionless factor to account for the altered speed of 
the planet at conjunction. Here, "+" refers to transits and 
"— " to occultations. One must also compute b using the 
eccentricity-dependent equations (|7]|8]i. 

For an eccentric orbit. Ting and Togi are generally unequal 
because the projected speed of the planet varies between 



ingress and egress. In practice the difference is slight; to 
leading order in R^,/a and e, 

Il^^ecosu;(^X{l-b'f\ (17) 



which is <10^2 g for a close-in planet with i?*/a = 0.2, 
and even smaller for more distant planets. For this reason 
we will use a single symbol r to represent either the ingress 
or egress duration. Another important timescale is T = 
Ttot — the interval between the halfway points of ingress 
and egress (sometimes referred to as contact times 1 .5 and 
3.5). 



In the limits e 



0, i?p < i?^ < a, and 6 < 1 - fc 



(which excludes near-grazing events), the results are greatly 
simplified: 



T^ToVl-b^, 



Tpk 



where Tn is the characteristic time scale 



Tn 



R.P 



13 hr 



P 

1 yr 



1/3 



^1/3 



(18) 



(19) 



For eccentric orbits, the additional factor given by equa- 
tion (fT6] l should be applied. Note that in deriving equa- 
tion ( fT9l ), we used Kepler's third law and the approximation 
Mp ^ AI^, to rewrite the expression in terms of the stel- 
lar mean density p*. This is a hint that eclipse observations 
give a direct measure of p*, a point that is made more ex- 
plicit in Section 3.1. 

2.4 Loss of light during eclipses 

The combined flux F{t) of a planet and star is plotted 
in Figure [T] During a transit, the flux drops because the 
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planet blocks a fraction of the starlight. Then the flux rises 
as the planet's dayside comes into view. The flux drops 
again when the planet is occulted by the star. Conceptually 
we may dissect F{t) as 



F{t) = F^{t)+Fp{ty 



Pau-e,{t)Fi,{t) transits, 







outside eclipses, 
Oiocc{t)Fp{t) occultations. 

(20) 

where Fi, , Fp are the fluxes from the stellar and planetary 
disks, and the a's are dimensionless functions of order unity 
depending on the overlap area between the stellar and plan- 
etary disks. In general F^, may vary in time due to flares, 
rotation of star spots and plages, rotation of the tidal bulge 
raised by the planet, or other reasons, but for simplicity of 
discussion we take it to be a constant. In that case, only the 
ratio f{t) = F{t)/ Fi, is of interest. If we let Ip and F be 
the disk-averaged intensities of the planet and star, respec- 
tively, then Fp/Fi, = k'^Ip/F and 



fit) = 1 + fc 



fc at 




it) 



transits, 

outside eclipses, 
occultations. 

(21) 

Time variations in Ip are caused by the changing illumi- 
nated fraction of the planetary disk (its phase function), as 
well as any changes intrinsic to the planetary atmosphere. 
To the extent that Ip is constant over the relatively short 
timespan of a single eclipse, all of the observed time vari- 
ation is from the a functions. As a starting approximation 
the a's are trapezoids, and f{t) is specified by the depth 6, 
duration T, ingress or egress duration r, and time of con- 
junction tc, as shown in Figure|2l For transits the maximum 
loss of light is 
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Ip{ttia) 

F 



(22) 



and in the usual case when the light from the planetary 
nightside is negligible, ^tia ~ fc^- For occultations. 



! ^pif occ) 

F ■ 



(23) 



In the trapezoidal approximation the flux variation dur- 
ing ingress and egress is linear in time. In reality this is 
not true, partly because of the nonuniform motion of the 
stellar and planetary disks. More importantly, even with 
uniform motion the overlap area between the disks is not a 
linear function of time [see equation (1) of Mandel & Agol 
(2002)]. In addition, the bottom of a transit light curve is 
not flat because real stellar disks do not have uniform inten- 
sity, as explained in the next section. 

2.5 Limb darkening 

Real stellar disks are brighter in the middle and fainter 
at the edge (the limb), a phenomenon known as limb dark- 
ening. This causes the flux decline during a transit to be 



larger than fc^ when the planet is near the center of the star, 
and smaUer than k^ when the planet is near the limb. The 
effect on the light curve is to round off the bottom and blur 
the 2nd and 3rd contact points, as shown in Figure H) Limb 
darkening is a consequence of variations in temperature and 
opacity with altitude in the stellar atmosphere. The sight- 
line to the limb follows a highly oblique path into the stel- 
lar atmosphere, and therefore an optical depth of unity is 
reached at a higher altitude, where the temperature is cooler 
and the radiation is less intense. The resulting intensity pro- 
file I{X, Y) is often described with a fitting formula such 
as 

/oc 1 -/i) -U2(l (24) 



where /i = y/1 — — and {^1,1/2} are constant co- 
efficients that may be calculated from stellar-atmosphere 
models or measured from a sufficiently precise transit light 
curve. The decision to use the quadratic function of equa- 
tion (l24l i or another of the various limb-darkening "laws" 
(better described as fitting formulas) is somewhat arbi- 
trary. Claret (2004) provides a compilation of theoretical 
coefficients, and advocates a four-parameter law. South- 
worth (2008) investigates the results of fitting the same data 
set with different limb-darkening laws. Mandel & Agol 
(2002) give accurate expressions for at^^{t) for some limb- 
darkening laws, and Gimenez (2006) shows how to com- 
pute atiii{t) for an arbitrary law based on earlier work by 
Kopal(1979). 

By using one of these methods to calculate the flux of 
a limb-darkened disk with a circular obstruction, it is usu- 
ally possible to model real transit light curves to within the 
measurement errors (see Section 4.3). In principle, calcu- 
lations of occultation light curves should take the planetary 
limb darkening into account, though the precision of cur- 
rent data has not justified this level of detail. Likewise, in 
exceptional cases it may be necessary to allow for depar- 
tures from circular shapes, due to rotational or tidal defor- 
mation (see Section 6). Modelers of eclipsing binary stars 
have long needed to take into account these and other subtle 
effects (Kallrath & Milone 2009, Hilditch 2001). 

More generally, the loss of light depends on the intensity 
of the particular patch of the photosphere that is hidden by 
the planet. The planet provides a raster scan of the stellar 
intensity across the transit chord. In this manner, star spots 
and plages can be detected through the flux anomalies that 
are observed when the planet covers them. (An example is 
given in the top right panel of Figure[8]) Even spots that are 
not along the transit chord can produce observable effects, 
by causing variations in and thereby causing Jtra to vary 
from transit to transit. 

3. SCIENCE FROM ECLIPSES 

3.1 Determining absolute dimensions 

Chief among the reasons to observe transits is to de- 
termine the mass and radius of the planet. Ideally one 
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Fig. 4. — Transits of tiie giant planet HD 209458b observed at 
wavelengths ranging from 0.32 /im (bottom) to 0.97 (top). At 
shorter wavelengths, the limb darkening of the star is more pro- 
nounced, and the bottom of the light curve is more rounded. The 
data were collected with the Hubble Space Telescope by Knutson 
et al. (2007a). 



would like to know the mass in kilograms, and the radius 
in kilometers, to allow for physical modeling and compar- 
isons with solar system planets. With only a transit light 
curve, this is impossible. The light curve by itself reveals 
the planet-to-star radius ratio k = Rp/R^ w ^/S but not 
the planetary radius, and says nothing about the planetary 
mass. 

To learn the planetary mass, in addition to the light curve 
one needs the radial-velocity orbit of the host star, and 
in particular the velocity semiamplitude i^T*. Using equa- 
tion (66) of the chapter by Murray and Correia, and Kepler's 
third law, we may write 



Sim 



P 
2^ 



1/3 



(25) 



The observation of transits ensures sini « 1, and thereby 
breaks the usual Mp sin i degeneracy. However, the plan- 
etary mass still cannot be determined independently of the 
stellar mass. In the usual limit Mp ^ A/*, the data deter- 



f2/3 



but not Mp itself. 



mine Mp/Mf 

To determine the absolute dimensions of the planet, one 
must supplement the transit photometry and radial-velocity 
orbit with some external information about the star De- 
pending on the star, the available information may include 
its luminosity, spectral type, and spectrally-derived photo- 
spheric properties (effective temperature, surface gravity, 
and metallicity). A typical approach is to seek consistency 
between those data and stellar-evolutionary models, which 
relate the observable properties to the stellar mass, radius, 
composition, and age. In special cases it may also be pos- 
sible to pin down the stellar properties using interferometry 
(Baines et al. 2009), asteroseismology (Stello et al. 2009), 
or an eclipsing companion star. 

Besides 5, the transit light curve offers the observables 
Ttot and Tfun (or T and t), which can be used to solve for 
the impact parameter h and the scaled stellar radius Ri,/a. 
For non-grazing transits, in the limit i?p ^ i?* ^ a we 
may invert equations (11411151 ) to obtain the approximate for- 
mulas 



62 



R* 
a 



1 - (Tf„n/Ttot)2 



2(51/4 



P 



e sm uj 



(26) 
(27) 



If in addition t <^ T, such as is the case for small planets 
on non-grazing trajectories, the results are simplified still 
further to 



1-^/6 



T 



R^ 

a 



T 

/TV 



(51/4 p 



1 + e sin UJ 



(28) 
(29) 



The orbital inclination i may then be obtained using equa- 
tion These approximations are useful for theoretical 
calculations and for developing an intuition about how the 
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system parameters affect the observable light curve. For 
example, i?*/a controls the product of T and r while b 
controls their ratio. However, as mentioned earlier, for fit- 
ting actual data one needs a realistic limb-darkened model, 
Unked to a Keplerian orbital model. 

The dimensionless ratios R^^/a and Rp/a are impor- 
tant for several reasons: (i) They set the scale of tidal in- 
teractions between the star and planet, (ii) Rp/a deter- 
mines what fraction of the stellar luminosity impinges on 
the planet, as discussed in Section 3.4. (iii) i?*/a can be 
used to determine a particular combination of the stellar 
mean density and planetary mean density pp: 



k^Pp = 



37r 
GP2 



a 

R^ 



(30) 



This can be derived from Kepler's third law (Seager & 
Mallen-Ornelas 2003). Since is usually small, the sec- 
ond term on the left side of equation ( [30l ) is often negligible 
and can be determined purely from transit photometry. 

This method for estimating pi, has proven to be a useful 
diagnostic in photometric transit surveys: a true transit sig- 
nal should yield a value of p^, that is consistent with expec- 
tations for a star of the given luminosity and spectral type. 
Furthermore, once a precise Ught curve is available, is a 
valuable additional constraint on the stellar properties. 

Interestingly, it is possible to derive the planetary surface 
gravity gp = GMp/Rp independently of the stellar proper- 
ties: 

^ 27r VT^R\ 

P {Rp/a)^ Sim- ^ ' 

This is derived from equation (66) of the chapter by Murray 
and Correia and Kepler's third law (Southworth et al. 2007). 

In short, precise transit photometry and Doppler ve- 
locimetry lead to correspondingly precise values of the stel- 
lar mean density p^, and planetary surface gravity gp. How- 
ever, the errors in Mp and Rp are ultimately limited by the 
uncertainties in the stellar properties. 

3.2 Timing of eclipses 

The orbital period P can be determined by timing a se- 
quence of transits, or a sequence of occultations, and fitting 
a linear function 



(32) 



where tc[n] is the time of conjunction of the nth event. The 
times must first be corrected to account for the Earth's or- 
bital motion and consequent variations in the light travel 
time. When comparing transit and occultation times, one 
must further correct for the light travel time across the line- 
of-sight dimension of the planetary orbit. As long as there is 
no ambiguity in n, the error in P varies inversely as the to- 
tal number of eclipses spanned by the observations, making 
it possible to achieve extraordinary precision in P. 

If the orbit does not follow a fixed ellipse — due to forces 
from additional bodies, tidal or rotational bulges, general 



relativity, or other non-Keplerian effects — then there will 
be variations in the interval between successive transits, as 
well as the interval between transits and occultations and 
the shape of the transit light curve. These variations may 
be gradual parameter changes due to precession (Miralda- 
Escude 2002), or short-term variations due to other planets 
(Holman & MuiTay 2005, Agol et al. 2005) or moons (Kip- 
ping 2009). The effects can be especially large for bodies 
in resonant orbits with the transiting planet. By monitoring 
transits one might hope to detect such bodies, as discussed 
in the chapter by Fabrycky. 

When transits and occultations are both seen, a power- 
ful constraint on the shape of the orbit is available. For a 
circular orbit, those events are separated in time by P/2, 
but more generally the time interval depends on e and uj. 
To first order in e, integrating dt/df between conjunctions 
gives 
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P 
2" 



1 



(33) 



In this case, the timing of transits and occultations gives an 
estimate of ecosa;. Likewise the relative durations of the 
transit and the occultation depend on the complementary 
parameter e sin oj, 
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e smu! 



(34) 



Sterne (1940) and de Kort (1954) give the lengthy exact re- 
sults for arbitrary e and i. Because the uncertainty in Atc/P 
is typically smaller than that in Ttm/Tocc (by a factor of 
P/T), the eclipse data constrain ecosw more powerfully 
than e sincj. 

The resulting bounds on e are often valuable. For exam- 
ple, planets on close-in eccentric orbits are internally heated 
by the friction that accompanies the time-variable tidal dis- 
tortion of the planet. Empirical constraints on e thereby 
help to understand the thermal structure of close-in planets. 
For more distant planets, bounds on e are helpful in the sta- 
tistical analysis of exoplanetary orbits. As described in the 
chapter by Gumming and in Part V of this volume, the ob- 
served eccentricity distribution of planetary orbits is a clue 
about the processes of planet formation and subsequent or- 
bital evolution. 

Eclipse-based measurements of <tra, tocc, and P are al- 
most always more precise than those based on spectroscopic 
or astrometric orbital data. The eclipse-based results can 
greatly enhance the analysis of those other data. For ex- 
ample, the usual radial-velocity curve has 6 parameters, but 
if ^tia, ^occ, and P are known from eclipses, the number of 
free parameters is effectively reduced to 3, thereby boosting 
the achievable precision in the other 3 parameters. 

3.3 Transmission spectroscopy 

We have been implicitly assuming that the planetary sil- 
houette has a sharp edge, but in reality the edge is fuzzy. For 
gas giant planets there is no well-defined surface, and even 
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planets with solid surfaces may have thick atmospheres. 
During a transit, a small portion of the starlight will be fil- 
tered through the upper atmosphere of the planet, where it is 
only partially absorbed. The absorption will be wavelength- 
dependent due to the scattering properties of atoms and 
molecules in the planetary atmosphere. At the wavelength 
of a strong atomic or molecular transition, the atmosphere is 
more opaque, and the planet's effective silhouette is larger. 
This raises the prospect of measuring the transmission spec- 
trum of the planet's upper atmosphere and thereby gaining 
knowledge of its composition. 

To calculate the expected signal one must follow the ra- 
diative transfer of the incident starlight along a grazing tra- 
jectory through the planet's stratified atmosphere. The cal- 
culation is rather complicated (see, e.g., Seager & Sasselov 
2000, Brown 2001) but the order of magnitude of the effect 
is easily appreciated. For a strong transition, the effective 
size of the planet grows by a few atmospheric scale heights 
H, where 



H = 



(35) 



and T is the temperature, is the mean molecular mass, 
g is the local gravitational acceleration, and fcs is Boltz- 
mann's constant. Defining Rp as the radius within which 
the planet is optically thick at all wavelengths, the extra ab- 
sorption due to the optically thin portion of the atmosphere 
causes the transit depth to increase by 



AS=^-^^^±^-^-S^2N,S[^], (36) 



nRl 



ttRI 



H 

R, 



where Nh, the number of scale heights, is of order unity. 
The signal is most readily detectable for planets with large 
H: low surface gravity, low mean molecular mass, and high 
temperature. For a "hot Jupiter" around a Sun-like star {5 = 
0.01, T « 1300 K, .9 « 25 m s^^, /i,„ = 2 amu) the signal 
is A(5 ^ 10^^. For an Earth-like planet around a Sun-like 



star {5 = 10-'*, T w 273 K, .g w 10 m s" 



28 amu). 



the signal is A(5 10"''. 

The signal can be detected by observing a transit light 
curve at multiple wavelengths, using different filters or a 
spectrograph. One then fits a limb-darkened light-curve 
model to the time series obtained at each wavelength, re- 
quiring agreement in the orbital parameters and allowing a 
value of 5 specific to each wavelength. The resulting vari- 
ations in (5(A) are expected to be of the order of magnitude 
given by equation ( |36] |. It is best to gather all of the data 
at the same time, because intrinsic stellar variability is also 
chromatic. 

3.4 Occultation spectroscopy 

As discussed in Section 2.4, when the planet is com- 
pletely hidden the starlight declines by a fraction (5occ = 
fc^Jp//*, where k is the planet-to-star radius ratio and /p//* 
is the ratio of disk-averaged intensities. Observations span- 
ning occultations thereby reveal the relative brightness of 



the planetary disk, if k is already known from transit obser- 
vations. The planetary radiation arises from two sources: 
thermal radiation and reflected starlight. Because the planet 
is colder than the star, the thermal component emerges at 
longer wavelengths than the reflected component. 

For the moment we suppose that the planet is of uni- 
form brightness, and that the observing wavelength is long 
enough for thermal emission to dominate. Approximating 
the planet and star as blackbody radiators. 



<5occ(A) 



fc2 Bx{Tp) 



where B\ (T) is the Planck function, 

2hc^ 1 



Xb ghc/{XkBT) _ I 



2kBT 



(37) 



, (38) 



in which T is the temperature, A is the wavelength, h is 
Planck's constant, and c is the speed of light. The limiting 
cases are for the "Rayleigh-Jeans" limit A ^ hc/{kBT). 
The decrement 5occ that is observed by a given instrument 
is obtained by integrating equation ( |37] | over the bandpass. 

Even when the planetary radiation is not described by the 
Planck law, one may define a brightness temperature Tb(A) 
as the equivalent blackbody temperature that would lead to 
the observed value of (5occ(A). The brightness temperature 
is sometimes a convenient way to describe the wavelength- 
dependent intensity even when it is not thermal in origin. 

There may be departures from a blackbody spectrum — 
spectral features — discernible in the variation of brightness 
temperature with wavelength. In contrast with transmission 
spectroscopy, which refers to starlight that grazes the plane- 
tary limb (terminator), here we are referring to the emission 
spectrum of the planet averaged over the visible disk of the 
dayside. Occultation spectroscopy and transit spectroscopy 
thereby provide different and complementary information 
about the planetary atmosphere. 

It is possible to measure the reflectance spectrum of the 
planet's dayside by observing at shorter wavelengths, or ac- 
curately subtracting the thermal emission. The occultation 
depth due to reflected light alone is 



<5occ(A) = Ax 



Rc 



(39) 



where Ax is the geometric albedo, defined as the flux re- 
flected by the planet when viewed at opposition (full phase), 
divided by the flux that would be reflected by a flat and per- 
fectly diffusing surface with the same cross-sectional area 
as the planet. One of the greatest uncertainties in atmo- 
spheric modeling is the existence, prevalence, and compo- 
sition of clouds. Since clouds can produce very large albedo 
variations, reflectance spectroscopy may help to understand 
the role of clouds in exoplanetary atmospheres. 

For a close-in giant planet, the reflectance signal is 
~10~^ while for an Earthlike planet at 1 AU it is ^10"^. 
The detection prospects are better for closer-in planets. 
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However, the closest-in planets are also the hottest, and 
their radiation may be dominated by thermal emission. An- 
other consideration is that planets with T > 1500 K are ex- 
pected to be so hot that all potentially cloud-forming con- 
densable materials are in gaseous form. The theoretically 
predicted albedos are very low, of order 10^"^, due to strong 
absorption by neutral sodium and potassium. 

Real planets do not have uniformly bright disks. Gaseous 
planets are limb darkened or brightened, and may have 
latitudinal zones with high contrast, like Jupiter. Rocky 
planets may have surface features and oceans. To the ex- 
tent that departures from uniform brightness could be de- 
tected, occultation data would provide information on the 
spatially-resolved planetary dayside. Specifically, the light 
curve of the ingress or egress of an occultation gives a 
one-dimensional cumulative brightness distribution of the 
planet. 

3.5 The Rossiter-McLaughlin effect 

In addition to the spectral variations induced by the plan- 
etary atmosphere, there are spectral variations arising from 
the spatial variation of the stellar spectrum across the stellar 
disk. The most pronounced of these effects is due to stel- 
lar rotation: light from the approaching half of the stellar 
disk is blueshifted, and light from the receding half is red- 
shifted. Outside of transits, rotation broadens the spectral 
lines but does not produce an overall Doppler shift in the 
disk-integrated starlight. However, when the planet covers 
part of the blueshifted half of the stellar disk, the integrated 
starlight appears slightly redshifted, and vice versa. 

Thus, the transit produces a time-variable spectral dis- 
tortion that is usually manifested as an "anomalous" radial 
velocity, i.e., a Doppler shift that is greater or smaller than 
the shift expected from only the star's orbital motion. Fig- 
ure |5] illustrates this effect. It is known as the Rossiter- 
McLaughlin (RM) effect, after the two astronomers who 
made the first definitive observations of this kind for binary 
stars, in 1924. 

The maximum amplitude of the anomalous radial veloc- 
ity is approximately 

AFrm ~ fcVl-^'K sin^,), (40) 

where Vi, sin i-^ is the line-of-sight component of the 
stellar equatorial rotation velocity. For a Sun-like star 
{Vi, sin ii,=2 km s~^), the maximum amplitude is ^20 m s^^ 
for a Jovian planet and ^0.2 m s^^ for a terrestrial planet. 
This amplitude may be comparable to (or even larger than) 
the amplitude of the spectroscopic orbit of the host star. 
Furthermore it is easier to maintain the stability of a spec- 
trograph over the single night of a transit than the longer 
duration of the orbital period. Hence the RM effect is an 
effective means of detecting and confirming transits, pro- 
viding an alternative to photometric detection. 

In addition, by monitoring the anomalous Doppler shift 
throughout a transit, it is possible to measure the angle on 



the sky between the planetary orbital axis and the stellar 
rotation axis. Figure |6] shows three trajectories of a transit- 
ing planet that have the same impact parameter, and hence 
produce identical light curves, but that have different ori- 
entations relative to the stellar spin axis, and hence pro- 
duce different RM signals. The signal for a well-aligned 
planet is antisymmetric about the midtransit time (left pan- 
els), whereas a strongly misaligned planet that blocks only 
the receding half of the star will produce only an anomalous 
blueshift (right panels). 

A limitation of this technique is that it is only sensitive 
to the angle between the sky projections of the spin and or- 
bital angular momentum vectors. The true angle between 
those vectors is usually poorly constrained because is un- 
known. Nevertheless it may be possible to tell whether the 
planetary orbit is prograde or retrograde, with respect to the 
direction of stellar rotation. It is also possible to combine 
results from different systems to gain statistical knowledge 
about spin-orbit alignment. 

More broadly, just as transit photometry provides a raster 
scan of the intensity of the stellar photosphere along the 
transit chord, RM data provide a raster scan of the line-of- 
sight velocity field of the photosphere. This gives an inde- 
pendent measure of the projected rotation rate u^, sin i^,, and 
reveals the velocity structure of starspots or other features 
that may exist on the photosphere. 

4. OBSERVING ECLIPSES 

4.1 Discovering eclipsing systems 

Figure|2]shows the rate of exoplanet discoveries over the 
last 20 years, highlighting the subset of planets known to 
transit. Transits have been discovered in two ways. One 
way is to find the planet in a Doppler survey (see chap- 
ter by Lovis and Fischer), and then check for transits by 
monitoring the brightness of the star throughout the inferior 
conjunction. This was the path to discovery for 6 of the 64 
known systems, including the first example, HD 209458b 
(Charbonneau et al. 2000, Henry et al. 2000, Mazeh et 
al. 2000). The probability that transits will occur is given 
by equation (9) and the times of inferior conjunction can be 
calculated from the parameters of the spectroscopic orbit 
(Kane 2007). 

The other way is to conduct photometric surveillance of 
stars that are not yet known to have any planets, an idea dat- 
ing back to Struve (1952). As an illustration let us consider 
the requirements for a program to discover hot Jupiters. As 
discussed in Section 2.2, a bare minimum of lO'^ Sun-like 
stars must be examined to find a transiting planet with an 
orbital distance of 0.05 AU. The observations must be pre- 
cise enough to detect a 1 % flux drop, and must extend for 
at least several times longer than the 3 day orbital period. 
The difficulty increases dramatically for smaller planets in 
wider orbits: in an idealized survey of nearby field stars 
limited by photon noise (see Section 4.2) the sensitivity to 
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Fig. 5. — Illustration of the Rossiter-McLaughlin (RM) effect. The three columns show three successive phases of a transit. The 
first row shows the stellar disk, with the grayscale representing the projected rotation velocity: the approaching limb is black and the 
receding limb is white. The second row shows the corresponding stellar absorption line profiles, assuming rotation to be the dominant 
broadening mechanism. The "bump" occurs because the planet hides a fraction of the light that contributes a particular velocity to the 
line-broadening kernel. The third row shows the case for which other line-broadening mechanisms are important; here the RM effect is 
manifested only as an "anomalous Doppler shift." Adapted from Gaudi & Winn (2007). 
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Fig. 6. — Using the RM effect to measure the angle A between the sky projections of the orbital and stellar-rotational axes. Three 
different possible trajectories of a transiting planet are shown, along with the corresponding RM signal. The trajectories all have the 
same impact parameter and produce the same light curve, but they differ in A and produce different RM curves. The dotted lines are for 
the case of no limb darkening, and the solid lines include limb darkening. From Gaudi & Winn (2007). 
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Fig. 7. — Rate of exoplanet discovery. The ligtit bars show the number of announcements of newly discovered planets, and the dark bars 
show the subset of transiting planets. Compiled with data from exoplanet.eu, an encyclopedic web site maintained by J. Schneider. 



transiting planets scales approximately as R^/P^/^ (Gaudi 
2005; Gaudi, Seager, & Mallen-Ornelas 2005). More real- 
istic calculations take into account the spread in radius and 
luminosity among the surveyed stars, the duration and time- 
sampling of observations, the appropriate threshold value of 
the signal-to-noise ratio, and other factors (see, e.g.. Pepper 
et al. 2003, Beatty & Gaudi 2008). 

Photometric surveys are much less efficient than Doppler 
surveys, in the sense that only a small fraction of plan- 
ets transit, and even when transits occur they are under- 
way only a small fraction of the time. On the other hand, 
the starting equipment for a photometric survey is mod- 
est, consisting in some cases of amateur-grade telescopes 
or telephoto lenses and cameras, whereas Doppler surveys 
require upfront a large telescope and sophisticated spectro- 
graph. For these reasons the royal road of eclipses enticed 
many astronomers to embark on photometric transit sur- 
veys. More than a dozen surveys were undertaken, includ- 
ing a few longitudinally-distributed networks to provide 
more continuous time coverage. Major efforts were made 
to automate the observations, and to develop algorithms for 
precise wide-field photometry (Tamuz et al. 2005, Kovacs 
et al. 2005) and transit detection (Kovacs et al. 2002). 
Home (2003) took stock of all the ongoing and planned 
surveys, and ventured to predict a bounty of 10-100 new 
planets per month, in an article subtitled "Hot Jupiters Ga- 
lore." 

The royal road turned out to have some potholes. One 
obstacle was the high rate of "false positives," signals that 
resemble planetary transits but are actually grazing eclipses 
of a binary star, or an unresolved combination of an ecfips- 
ing binary star and a third star. In the latter case, the deep 
eclipses of the binary star are diluted to planet-Uke propor- 
tions by the constant light of the third star. In some surveys 



the false positives outnumbered the planets by 10 to 1. Rul- 
ing them out required spectroscopy with large telescopes, 
which became the bottleneck in the discovery process (see, 
e.g., O'Donovan et al. 2007). Another obstacle was corre- 
lated noise ("red noise") in the survey photometry (see Sec- 
tion 4.3). Pont et al. (2006) showed that red noise slashed 
the sensitivity of the search algorithms, thereby providing a 
quantitative solution to the "Horne problem" of why tran- 
siting planets were not being found as rapidly as expected. 

Only a few of the surveyors were able to overcome these 
obstacles. The 5 most successful surveys were OGLE, 
which used aim telescope to survey 14- 16th magnitude 
stars; and the TrES, XO, HAT, and SuperWASP surveys, 
which used ssO.l m lenses to survey 10-12th magnitude 
stars. Many of today's transiting planets are named after 
these surveys: for more details see Udalski et al. (2002), 
Alonso et al. (2004), McCullough et al. (2005), Bakos et 
al. (2007), and PoUacco et al. (2006). The OGLE planets 
were found first, although the brighter stars from the wider- 
field surveys were far more amenable to false-positive re- 
jection and detailed characterization. 

On the other end of the cost spectrum are photometric 
surveys performed from space. By avoiding the deleterious 
effects of the Earth's atmosphere, it is possible to beat the 
precision of ground-based observations (see Section 4.2). 
It is also possible to avoid the usual interruptions due to 
the vagaries of the weather, and even to avoid the day-night 
cycle if the spacecraft is in an orbit far from Earth. The 
two ongoing space -based missions CoRoT and Kepler are 
described in Sections 5.1 and 6, respectively. 

One space-based survey that did not discover any plan- 
ets neverthless produced an interesting result. Gilliland et 
al. (2000) used the Hubble Space Telescope to seek close-in 
giant planets in the globular cluster 47 Tucanae. No transits 
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were found, even though 17 were expected based on the 
survey characteristics and the observed frequency of close- 
in giant planets around nearby field stars. The absence of 
close-in planets could be due to the crowded stellar envi- 
ronment (which could inhibit planet formation or disrupt 
planetary systems) or the cluster's low metallicity (which 
has been found to be correlated with fewer planets in the 
Doppler surveys; see the chapters by Lovis and Fischer 
and by Gumming). The null results of a subsequent sur- 
vey of the less-crowded outskirts of 47 Tucanae suggest that 
crowding is not the primary issue (Weldrake et al. 2005). 

4.2 Measuring the photometric signal 

Even after a transiting planet is discovered, it takes a 
careful hand to measure the transit signal precisely enough 
to achieve the scientific goals set forth in Section 3. The loss 
of light is only 1 % for a Sun-like star crossed by a Jupiter- 
sized planet, and 0.01% for an Earth-sized planet. Occul- 
tations produce still smaller signals. This is the domain of 
precise time-series differential photometry. The term "dif- 
ferential" applies because only the fractional variations are 
of interest, as opposed to the actual intensity in Janskys or 
other standardized units. Eclipses can also be detected via 
spectroscopy (Section 3.5) or polarimetry (Carciofi & Ma- 
gaUiaes 2005) but here we focus on the photometric signal, 
with emphasis on ground-based optical observations. 

First you must know when to observe. To observe an 
eclipse requires a triple coincidence: the eclipse must be 
happening, the star must be above the horizon, and the Sun 
must be down. Transit times can be predicted based on a 
sequence of previously measured transit times, by fitting 
and extrapolating a straight line (equation l32Ti. Occultation 
times can also be predicted from a listing of transit times, 
but are subject to additional uncertainty due to the depen- 
dence on e and lo (equation[33]l. 

Next you should monitor the flux of the target star along 
with other nearby stars of comparable brightness. The mea- 
sured fluxes are affected by short-term variations in atmo- 
spheric transparency, as well as the gradual change in the 
effective atmospheric path length (the airmass) as the star 
rises and sets. However, the ratios of the fluxes between 
nearby stars are less affected. As long as some of the com- 
parison stars are constant in brightness then the relative flux 
of the target star can be tracked precisely. 

This task is usually accomplished with a charge-coupled 
device (CCD) imaging camera and software for calibrat- 
ing the images, estimating the sky background level, and 
counting the photons received from each star in excess of 
the sky level {aperture photometry). Howell (1999) ex- 
plains the basic principles of CCDs and aperture photom- 
etry. For more details on differential aperture photome- 
try using an ensemble of comparison stars, see Gilliland & 
Brown (1988), Kjeldsen & Frandsen (1992), and Everett & 
Howell (2001). In short, the fluxes of the comparison stars 
are combined, and the flux of the target star is divided by 
the comparison signal, giving a time series of relative flux 



measurements spanning the eclipse with as little noise as 
possible. For bright stars observed at optical wavelengths, 
among the important noise sources are photon noise, scintil- 
lation noise, differential extinction, and flat-fielding errors, 
which we now discuss in turn. 

Photon noise refers to the unavoidable fluctuations in the 
signal due to the quantization of light. It is also called Pois- 
son noise because the photon count rate obeys a Poisson 
distribution. If a star delivers N photons s^^ on average, 
then the standard deviation in the relative flux due to photon 
noise is approximately (iVAf)~^/^ for an exposure lasting 
At seconds. This noise source affects the target and com- 
parison stars independently. The sky background also in- 
troduces Poisson noise, which can be troublesome for faint 
stars or infrared wavelengths. The photon noise in the com- 
parison signal can be reduced by using many bright com- 
parison stars. Beyond that, improvement is possible only by 
collecting more photons, using a bigger telescope, a more 
efficient detector, or a wider bandpass. 

Scintillation is caused by fluctuations in the index of re- 
fraction of air The more familiar term is "twinkling." For 
integration times At ^ 1 s, the standard deviation in the rel- 
ative flux due to scintillation is expected to scale with tele- 
scope diameter D, observatory altitude h, integration time, 
and airmass as 

(Airmass) / h \ 

-scin = -o^^575(A^ {-^^) ' (41) 

based on a theory of atmospheric turbulence by Reiger (1963), 
with empirical support from Young (1967) and others. The 
coefficient ctq is often taken to be 0.064 when D is ex- 
pressed in centimeters and At in seconds, but this must be 
understood to be approximate and dependent on the local 
meteorology. Scintillation affects both the target and com- 
parison stars, although for closely-spaced stars the vari- 
ations are correlated (Ryan & Sandler 1998). One also 
expects scintillation noise to decrease with wavelength. 
Thus, scintillation noise is reduced by employing a large 
telescope (or combining results from multiple small tele- 
scopes), choosing nearby comparison stars, and observing 
at a long wavelength from a good site. 

Differential extinction is used here as a shorthand for 
"second-order color-dependent differential extinction." To 
first order, if two stars are observed simultaneously, their 
fluxes are attenuated by the Earth's atmosphere by the same 
factor and the flux ratio is preserved. To second order, the 
bluer star is attenuated more, because scattering and ab- 
sorption are more important at shorter wavelengths. This 
effect causes the flux ratio to vary with airmass as well as 
with short-term transparency fluctuations. It can be reduced 
by choosing comparison stars bracketing the target star in 
color, or using a narrow bandpass. The advantage of a nar- 
row bandpass must be weighed against the increased photon 
noise. 

Flat fielding is the attempt to correct for nonuniform illu- 
mination of the detector and pixel-to-pixel sensitivity varia- 
tions, usually by dividing the images by a calibration image 
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of a uniformly-lit field. If the stars were kept on the same 
pixels throughout an observation, then flat-fielding errors 
would not affect the flux ratios. In reality the light from a 
given star is detected on different pixels at different times, 
due to pointing errors, focus variations, and seeing varia- 
tions. Imperfect flat fielding coupled with these variations 
produce noise in the light curve. The impact of flat-fielding 
errors is reduced by ensuring the calibration images have 
negligible photon noise, maintaining a consistent pointing, 
and defocusing the telescope. Defocusing averages down 
the interpixel variations, reduces the impact of seeing vari- 
ations, and allows for longer exposure times without sat- 
uration, thereby increasing the fraction of the time spent 
collecting photons as opposed to resetting the detector. De- 
focusing is good for what ails you, as long as the stars do 
not blend together 

This discussion of noise sources is not exhaustive; it is in 
the nature of noise that no such listing can be complete. For 
example, the gain of the detector may drift with temper- 
ature, or scattered moonlight may complicate background 
subtraction. A general principle is to strive to keep every- 
thing about the equipment and the images as consistent as 
possible. Another good practice is to spend at least as much 
time observing the star before and after the eclipse as during 
the eclipse, to establish the baseline signal and to character- 
ize the noise. 

It is often advisable to use a long-wavelength bandpass, 
not only to minimize scintillation and differential extinc- 
tion, but also to reduce the effects of stellar limb darken- 
ing on the transit light curve. The degree of limb darken- 
ing diminishes with wavelength because the ratio between 
blackbody spectra of different temperatures is a decreas- 
ing function of wavelength. Transit light curves observed 
at longer wavelengths are "boxier," with sharper corners 
and flatter bottoms. All other things being equal, this re- 
duces the statistical uncertainties in the transit parameters, 
but other factors should also be considered. For example, at 
infrared wavelengths, limb darkening may be small, but the 
sky background is bright and variable. 

The space-based observer need not worry about scintil- 
lation and differential extinction, and enjoys a low back- 
ground level even at infrared wavelengths. Extremely pre- 
cise photometry is possible, limited only by the size of 
the telescope and the degree to which the detector is well- 
calibrated. With the Hubble Space Telescope a precision of 
approximately 10^"' per minute-long integration has been 
achieved, a few times better than the best ground-based 
light curves. Figure |8] allows for some side-by-side com- 
parisons of ground-based and space-based data. However, 
going to space is not a panacea. When unforeseen calibra- 
tion issues arise after launch, they are difficult to resolve. In 
low Earth orbit, there are also unhelpful interruptions due 
to occultations by the Earth, as well as problems with scat- 
tered Earthshine and cosmic rays. More remote orbits offer 
superior observing conditions, but require greater effort and 
expense to reach. 
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Fig. 8. — Examples of transit light curves based on ground-based 
observations (left) and space-based observations (right). In the top 
right panel, the "bump" observed just before midtransit is inter- 
preted as the covering of a dark starspot by the planet. From upper 
left to lower right, the references are Winn et al. (2007a), Rabus 
et al. (2009), Winn et al. (2007b), Pont et al. (2008), Holman et 
al. in prep., Gillon et al. (2007), Winn et al. (2008), and Carter et 
al. (2009). 
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4.3 Interpreting the photometric signal 

Once you have an eclipse light curve, the task remains 
to derive the basic parameters {S, T, r, tc} and their uncer- 
tainties, as well as the results for other quantities such as gp 
and p^, that can be derived from those parameters. The an- 
alytic equations given in Section 2 for the eclipse duration, 
timing, and other properties are rarely used to analyze data. 
Rather, a parametric model is fitted to the data, based on 
the numerical integration of Kepler's equation to calculate 
the relative positions of the star and planet, as well as one 
of the prescriptions mentioned in Section 2.4 for computing 
the loss of light from the limb-darkened stellar photosphere. 

The first task is writing a code that calculates the light 
curve of the star-planet system as a function of the orbital 
parameters and eclipse parameters. Then, this code is used 
in conjunction with one of many standard routines to op- 
timize the parameter values, typically by minimizing the 
sum-of-squares statistic. 



N r 



X 



E 



/,(obs) - /i(calc) 



(42) 



where /i(obs) is the observed value of relative flux at 
time ti, /i(calc) is the calculated flux (depending on the 
model parameters), and <Ti is the measurement uncertainty. 
The techniques found in standard works such as Numerical 
Recipes (Press et al. 2007) are applicable here, although a 
few points of elaboration are warranted for the specific con- 
text of eclipse photometry. 

It is common to adopt a Bayesian attitude, in which the 
parameters are viewed as random variables whose probabil- 
ity distributions ("posteriors") are constrained by the data. 
This can be done in a convenient and elegant fashion us- 
ing the Monte Carlo Markov Chain (MCMC) method, in 
which a chain of points is created in parameter space using 
a few simple rules that ensure the collection of points will 
converge toward the desired posterior This method gives 
the full multidimensional joint probability distribution for 
all the parameters, rather than merely giving individual er- 
ror bars, making it easy to visualize any fitting degeneracies 
and to compute posteriors for any combination of parame- 
ters that may be of interest. Although a complete MCMC 
briefing is beyond the scope of this chapter, the interested 
reader should consult the textbook by Gregory (2005) as 
well as case studies such as Holman et al. (2006), Collier 
Cameron et al. (2007), and Burke et al. (2007). 

A vexing problem is the presence of correlated noise, 
as mentioned in Section 4.2. The use of equation (|42] | is 
based on the premise that the measurement errors are sta- 
tistically independent. In many cases this is plainly false. 
Real light curves have bumps, wiggles, and slopes span- 
ning many data points. These can be attributed to differen- 
tial extinction, flat-fielding errors, or astrophysical effects 
such as starspots. Thus the number of truly independent 
samples is smaller than the number of data points, and the 
power to constrain model parameters is correspondingly re- 



duced. Ignoring the correlations leads to false precision in 
the derived parameters, but accounting for correlations is 
not straightforward and can be computationally intensive. 
Some suggestions are given by Pont et al. (2006) and Carter 
& Winn (2009a). 

The treatment of stellar limb darkening presents another 
unwelcome opportunity to underestimate the parameter un- 
certainties. It is tempting to adopt one of the standard limb- 
darkening laws (see Section 2.5) and hold the coefficients 
fixed at values deemed appropriate for the host star, based 
on stellar-atmosphere models. However, with precise data 
it is preferable to fit for the coefficients, or at least to allow 
for some uncertainty in the atmospheric models. In those 
few cases where the data have been precise enough to test 
the models, the models have missed the mark (Claret 2009). 

Although fitting data is a job for a computer, it is never- 
theless useful to have analytic formulas for the achievable 
precision in the eclipse parameters. The formulas are handy 
for planning observations, and for order-of-magnitude esti- 
mates of the observability (or not) of effects such as vari- 
ations in transit times and durations. The formulas given 
here are based on the assumptions that the data have uni- 
form time sampling At, and independent Gaussian errors a 
in the relative flux. A useful figure of merit is Q = \/NS/(t, 
where S is the transit depth and N is the number of data 
points obtained during the transit. A Fisher information 
analysis (essentially a glorified error propagation) leads to 
estimates for the la uncertainties in the transit parameters 
(Cai-ter et al. 2008): 



Or 



Q-^T^t/2T, 



(43) 
(44) 
(45) 
(46) 



which are valid when 6-^1, limb darkening is weak, and 
the out-of-transit flux is known precisely. In this case, at^ < 
ax < Or- Correlated errors and limb darkening cause these 
formulas to be underestimates. 



5. SUMMARY OF RECENT ACHIEVEMENTS 
5.1 Discoveries of transiting planets 

Discoveries of exoplanets, and transiting exoplanets in 
particular, have abounded in recent years (Figure IT). As 
of December 2009, approximately 64 transiting planets are 
known, representing 15% of the total number of exoplanets 
discovered. Figure |9] shows their masses, radii, and orbital 
periods. It is important to remember that these planets are 
not a randomly selected subset of exoplanets. The proper- 
ties of the ensemble have been shaped by powerful selection 
effects in the surveys that led to their discovery, favoring 
large planets in short-period orbits. 

Despite these selection effects, the known transiting 
planets exhibit a striking diversity. They span three or- 
ders of magnitude in mass, and one order of magnitude in 
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Orbital Period [days] Orbital Period [days] 

Fig. 9. — Masses, radii, and orbital periods of the transiting exoplanets. Upper left: Radius versus mass, on logarithmic scales. The 
dotted lines are loci of constant mean density (0.5, 1, and 2 g cm^'^, from top to bottom). Upper right: Same, but on linear scales, and 
with the axes restricted to highlight the gas giants that dominate the sample. Lower left: Mass versus orbital period, on a logarithmic 
scale. The two long-period outliers are HD 17156b (P = 21 d) and HD 80606b (P = 111 d). Lower right: Same, but on a linear scale, 
and with axes restricted to highlight the gas giants. The anticorrelation between mass and orbital period is evident. 
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radius. Most are gas giants, comparable in mass and radius 
to Jupiter. There are also two planets with sizes more like 
Neptune and Uranus, as well as two even smaller planets 
with sizes only a few times larger than Earth, a category 
that has come to be known as "super-Earths." 

The two transiting super-Earths (Figure [TOll were found 
with completely different strategies. The CoRoT (COnvec- 
tion. Rotation, and planetary Transits) team uses a satellite 
equipped with a 0.27 m telescope and CCD cameras to ex- 
amine fields of ~ 10,000 stars for a few months at a time, 
seeking relatively short-period planets (P < 20 d). Along 
with several giant planets they have found a planet of radius 
1 .7 i?0 in a 20-hour orbit around a G dwarf star, producing 
a transit depth of only 3.4 x 10^"* (Leger et al. 2009). This 
is smaller than the detection threshold of any of the ground- 
based surveys, demonstrating the advantage of space-based 
photometry. However, it proved difficult to spectroscopi- 
cally confirm that the signal is indeed due to a planet, be- 
cause the host star is relatively faint and chromospherically 
active (Queloz et al. 2009). 

Another project, called MEarth (and pronounced mirth), 
seeks transits of small planets from the ground by focusing 
on very small stars (M dwarfs), for which even a super- 
Earth would produce a transit depth of order 1%. Because 
such stars are intrinsically faint, one must search the whole 
sky to find examples bright enough for follow-up work. 
That is why MEarth abandoned the usual survey concept in 
which many stars are monitored within a single telescope's 
field of view. Instead they monitor M dwarf stars one at 
a time, using several 0.4 m telescopes. Using this strategy 
they found a planet of radius 2.7 i?© in a 1.6-day period 
around a star of radius 0.21 Rq (Charbonneau et al. 2009). 
The large transit depth of 1.3% invites follow-up observa- 
tions to study the planet's atmosphere. 

At the other extreme, several transiting objects have 
masses greater than 10 Af jup, reviving the old debate about 
what should and should not be considered a planet. The 
radii of these massive objects are not much larger than 
Jupiter's radius, in agreement with predictions that be- 
tween about 1-50 Af Jup the pressure due to Coulomb forces 
(which would give R cx M^/"^) and electron degener- 
acy pressure (i? cx M^^/'^) conspire to mute the mass- 
dependence of the radius. 

Almost all of the transiting planets have short orbital pe- 
riods (<10 days), due to the decline in transit probability 
with orbital distance (equation [T2]). Two conspicuous ex- 
ceptions are HD 17156b with P = 21 days (Barbieri et 
al. 2009), and HD 80606b with P = 111 days (Moutou 
et al. 2009, Garcia-Melendo & McCullough 2009, Fossey 
et al. 2009). Both of those planets were discovered by the 
Doppler method and found to transit through photometric 
follow-up observations. They were recognized as high- 
priority targets because both systems have highly eccentric 
orbits oriented in such a way as to enhance the probability 
of eclipses. Thus despite their long orbital periods, their pe- 
riastron distances are small (<0.1 AU) along with all of the 
other known transiting planets. 



5.2 Follow-up photometry and absolute dimensions 

Follow-up observations of transits have allowed the ba- 
sic transit parameters {5, T, r} to be determined to within 
1 % or better, and absolute dimensions of planets to within 
about 5%. Compilations of system parameters are given by 
Ton-es, Winn, & Holman (2008) and Southworth (2009). 
The orbital periods are known with 8 significant digits 
in some cases. Ground-based observations have achieved 
a photometric precision of 250 ppm per 1 min sample, 
through the techniques described in Section 4.2. 

A pathbreaking achievement was the Hubble Space Tele- 
scope light curve of HD 209458 by Brown et al. (2001), 
with a time sampling of 80 s and a precision of 1 10 ppm, 
without the need for comparison stars. With the Spitzer 
Space Telescope the photon noise is generally higher, but 
there are compensatory advantages: there is little limb- 
darkening at at mid-infrared wavelengths, and uninter- 
rupted views of entire events are possible because the satel- 
lite is not in a low-earth orbit. Among the most spectacular 
data yet obtained in this field was the 33 hr observation by 
Knutson et al. (2007b) of the K star HD 189733, spanning 
both a transit and occultation of its giant planet. The data 
were gathered with Spitzer at a wavelength of 8 /im, and 
are shown in Figure [TT] They not only provided extremely 
precise light curves of the transit and occultation, but also 
showed that the combined flux rises gradually in between 
those two events, demonstrating that the dayside is hotter 
than the nightside. 

Among the properties of the close-in giant planets, a few 
patterns have emerged. The planetary mass is inversely re- 
lated to the orbital period (Mazeh, Zucker, & Pont 2005). 
This anticorrelation has been variously attributed to selec- 
tion effects, tidal interactions, and thermal evaporation of 
the planetary atmosphere. There is also a positive correla- 
tion between the metallicity of the host star and the inferred 
"core mass" of the planet (Guillot et al. 2006), by which 
is meant the mass of heavy elements required in models of 
the planetary interior that agree with the observed mass and 
radius. It is tempting to interpret this latter correlation as 
support for the core-accretion theory of planet formation. 

A persistent theme in this field, and a source of contro- 
versy and speculation, is that several of the transiting giant 
planets have radii that are 10-50% larger than expected from 
models of hydrogen-helium planets, even after accounting 
for the intense stellar heating and selection effects (see the 
chapter by Correia and Laskar). Among the possible ex- 
planations for these "bloated" planets are tidal heating (see, 
e.g.. Miller et al. 2009); unknown atmospheric constituents 
that efficiently trap internal heat (Burrows et al. 2007); en- 
hanced downward convection allowing the incident stellar 
radiation to reach significant depth (Guillot & Showman 
2002); and inhibition of convection within the planet that 
traps internal heat (Chabrier & Baraffe 2007). 

Likewise, a few planets are observed to be smaller 
than expected for a hydrogen-helium planet with the ob- 
served mass and degree of irradiation. Some examples are 
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Fig. 10. — Transits of two different super-Earths, GJ 1214b (left) and CoRoT-7b (right). The planets are approximately the same size, 
but because GJ 1214b orbits a small star (spectral type M4.5V) its transit depth is much larger than that of CoRoT-7b, which orbits a 
larger star (G9V). References: Charbonneau et al. (2009), Leger et al. (2009). 
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Fig. 11. — The combined 8 /im brightness of the K star HD 189733 and its giant planet, over a 33 hr interval including a transit and an 
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HD 149026b (Sato et al. 2005) and HAT-P-3b (Torres et 
al. 2007). The favored interpretation is that these planets 
are enriched in elements heavier than hydrogen and helium, 
and the increased mean molecular weight leads to a larger 
overall density. Some degree of enrichment might be ex- 
pected because Jupiter and Saturn are themselves enriched 
in hevay elements relative to the Sun. However the two 
planets just mentioned would need to be enriched still fur- 
ther; in the case of HD 149026b, theoretical models suggest 
that the planet is 65% heavy elements by mass. 

5.3 Atmospheric physics 

Atmospheric spectra have been obtained for several gas 
giant planets, especially the two "hot Jupiters" with the 
brightest host stars, HD 209458b and HD 189733b. Both 
transmission and emission spectroscopy have been under- 
taken, mainly with the space telescopes Hubble and Spitzer. 
Figures 12 and 13 show examples of both transmission 
(transit) and emission (occupation) spectra. These obser- 
vations have stimulated much theoretical work on exoplan- 
etary atmospheres (see the chapters by Chambers and by 
D' Angelo et al. for more details than are presented here). 

This enterprise began with the optical detection of neu- 
tral sodium in a transmission spectrum of HD 209458b by 
Charbonneau et al. (2002). For that same planet, Vidal- 
Madjar et al. (2003) measured an enormous transit depth 
of (15 ± 4)% within an ultraviolet bandpass bracketing 
the wavelength of the Lyman-a transition, which they at- 
tributed to neutral hydrogen gas being blown off the planet. 
For HD 189733b, a rise in transit depth toward shorter 
wavelengths was observed by Pont et al. (2007) and inter- 
preted as Rayleigh scattering in the planet's upper atmo- 
sphere. 

Recent attention has turned to infrared wavelengths, 
where molecules make their imprint. In the emission spec- 
tra, departures from a blackbody spectrum have been inter- 
preted as arising from water, methane, carbon monoxide, 
and carbon dioxide (see, e.g., Grillmair et al. 2008, Swain 
et al. 2009). Absorption features in transmission spectra 
have been interpreted as arising from water and methane 
(see, e.g.. Swain et al. 2008, Tinetti et al. 2007). 

Despite these impressive achievements, many issues re- 
garding the atmospheres of hot Jupiters remain unsettled. 
Controversies arise because spectral features are observed 
with a low signal-to-noise ratio and are not always repro- 
ducible, as illustrated in Figure [12] This is understand- 
able, as the observers have pushed the instruments beyond 
their design specifications to make these demanding mea- 
surements. In addition, the theoretical interpretation of the 
spectra is in a primitive state. A published spectrum is 
usually accompanied by a model that fits the data, but left 
unanswered are whether the model is unique and how well- 
constrained are its parameters. Recent work by Madhusud- 
han & Seager (2009) addresses this problem. 

One theme that has arisen in the last few years is that 
some hot Jupiters have an inversion layer in their upper at- 



mospheres, within which the temperature rises with height 
instead of the usual decline. The evidence for an inver- 
sion layer is emission in excess of a blackbody spectrum 
between 4-8 /im, where excess absorption was expected 
due to water vapor The interpretation is that water is seen 
in emission because it exists in a hot, tenuous stratosphere. 

Hot stratospheres develop when starlight is strongly ab- 
sorbed by some species at low pressure (at high altitude) 
where the atmosphere does not radiate efficiently. The 
identity of this absorber in hot Jupiter atmospheres has 
been a topic of debate. Gaseous titanium oxide and vana- 
dium oxide are candidates (Hubeny et al. 2003), as are 
photochemically-produced sulfur compounds (Zahnle et 
al. 2009). Meanwhile observers are searching for corre- 
lations between the presence or absence of an inversion 
layer, the degree of stellar irradiation, the magnitude of 
the day-night temperature difference, and other observable 
properties. 

A few of the known planets have highly eccentric or- 
bits and small periastron distances, and are therefore sub- 
ject to highly variable stellar irradiation. Observers have 
monitored the planetary thermal emission following perias- 
tron passages, to gauge the amplitude and timescale of the 
thermal response to stellar heating. In the most extreme 
case, Laughlin et al. (2009) watched the effective tempera- 
ture of the giant planet HD 80606b rise from about 800 K 
to 1500 K after periastron passage, and inferred that the 
characteristic radiative timescale of the upper atmosphere 
is about 4.5 hr. 

Although the thermal emission from hot Jupiters has 
been detected in many cases, at infrared and optical wave- 
lengths, there has been no unambiguous detection of 
starlight reflected from the planetary atmosphere. The best 
resulting upper limit on the visual planetary albedo is 0.17 
(with 3(7 confidence), for the case of HD 209458b (Rowe 
et al. 2008). This rules out highly reflective clouds of the 
sort that give Jupiter its visual albedo of 0.5. However, the 
limits have not been of great interest to the theorists, who 
predicted all along that the visual albedos of hot Jupiters 
would be very small. 

5.4 Tidal evolution and migration 

Almost all of the known transiting planets are close 
enough to their parent stars for tidal effects to be impor- 
tant. The tidal bulges on the star and planet provide "han- 
dles" for the bodies to torque each other, and tidal friction 
slowly drains energy from the system (see the chapter by 
Fabrycky). For a typical hot Jupiter, the sequence of events 
is expected to be as follows: (i) Over ^10® yr, the planet's 
rotational period is synchronized with its orbital period, and 
its obliquity (the angle between its rotational and orbital 
angular momentum vectors) is driven to zero, (ii) Over 
^10^ yr, the orbit is circularized, (iii) After '-^lO^^ yr (i.e. 
not yet), the stellar rotational period is synchronized with 
the orbital period, the stellar obliquity is driven to zero, and 
the orbit decays, leading to the engulfment of the planet. 
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Fig. 12. — Transmission (transit) spectroscopy of tiie gas giant HD 189733b, using the Hubble Space Telescope. The symbols with 
errors bars are measurements of the effective planet-to-star radius ratio as a function of wavelength. The dip at 1.6 /im was interpreted 
as evidence for water, and the rise at 2.1 p,m. as evidence for methane (Swain et al. 2008). However, subsequent observations at 1.7 
and 1.9 /im, shown with darker symbols and gray bands, disagree with the earlier results and are consistent with a Rayleigh scattering 
model (solid and dashed curves). From Sing et al. (2009). 
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Fig. 13. — Occultation spectroscopy of the gas giant HD 189733b, using the Spitzer Space Telescope. The points show measurements of 
the flux density ratio of the planet and star as a function of wavelength. The smaller and finer-sampled points are based on observations 
with a dispersive spectrograph while the larger points are based on broadband filter photometry. The three lines show the outputs of 
model atmospheres with varying choices for the parameters P„, specifying the efficiency of heat transfer from day side to nightside, and 
Ke, specifying the opacity of a putative high-altitude absorbing species. From Grillmair et al. (2008). 
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The timescales for these processes are highly uncertain, and 
more complex histories are possible if there are additional 
planets in the system or if the planet's internal structure is 
strongly affected by tidal heating. 

Tidal circularization is implicated by the fact that plan- 
ets with orbital periods shorter than ^10 days tend to 
have smaller orbital eccentricities than longer-period plan- 
ets. This fact was akeady known from Doppler surveys, 
but for eclipsing planets the eccentricities can be measured 
more precisely (see Section 3.2). It is possible that the 
"bloating" of some of the close-in giant planets (Section 
5.3) is related to the heat that was produced during the cir- 
cularization process, or that may still be ongoing. 

As for tidal synchronization, observations with Spitzer 
have revealed planets with a cold side facing away from the 
star and a hot side facing the star (see FigurefTTli. This could 
be interpreted as evidence for synchronization, although it 
is not definitive, because it is also possible that the heat from 
the star is re-radiated too quickly for advection to homoge- 
nize the upper atmosphere of the planet. 

Tidal decay of the orbit would be observable as a grad- 
ual decline in the orbital period (Sasselov 2003). For most 
systems the theoretical timescale for tidal decay is much 
longer than the age of the star, and indeed no evidence for 
this process has been found. However in at least one case 
(WASP- 18b; Hellier et al. 2009) the theoretical timescale 
for tidal decay is much shorter than the stellar age, because 
of the large planetary mass and short orbital period. The 
existence of this system suggests that the theoretical expec- 
tations were wrong and dissipation is slower in reality. 

Likewise, in most cases one would not expect that 
enough time has elapsed for tides to modify the star's spin 
rate or orientation (Barker & Ogilvie 2009). This suggests 
that the measurements of the projected spin-orbit angle us- 
ing the RM effect (Section 3.5) should be interpreted in the 
context of planet formation and evolution rather than tides. 
A close spin-orbit alignment is expected because a star and 
its planets inherit their angular momentum from a common 
source: the protostellar disk. However, for hot Jupiters 
there is the complication that they presumably formed at 
larger distances and "migrated" inward through processes 
that are poorly understood. Some of the migration theo- 
ries predict that the original spin-orbit alignment should 
be preserved, while others predict occasionally large mis- 
alignments. For example, tidal interactions with the pro- 
toplanetary disk (see chapter by Showman et al.) should 
drive the system into close alignment (Marzari & Nelson 
2009), while planet-planet scattering or Kozai oscillations 
with tidal friction (see chapter by Fabrycky) should result 
in misaligned systems. 

The projected spin-orbit angle has been measured for 
about 20 exoplanets, all of them close-in giants. Some ex- 
amples of data are shown in Figure [14] In many cases the 
results are consistent with good alignment, with measure- 
ment precisions ranging from 1-20° (for a recent summary, 
see Fabrycky & Winn 2009). However there are now at 
least 4 clear cases of misaligned systems. One such case is 



XO-3b, a massive planet in a close-in eccentric orbit that is 
tilted by more than 30° with respect to the stellar equator 
(Hebrard et al. 2008, Winn et al. 2009a). Even more dra- 
matic is HAT-P-7b, for which the planetary orbit and stellar 
spin axis are tilted by more than 86° (Winn et al. 2009b, 
Narita et al. 2009). The planetary orbit is either polar (go- 
ing over the north and south poles of the star) or retrograde 
(revolving in the opposite direction as the star is rotating). 
Another possible retrograde system is WASP- 17b (Ander- 
son et al. 2009). 

6. FUTURE PROSPECTS 

Eclipses are the "here and now" of exoplanetary sci- 
ence. It seems that every month, someone reports a startling 
observational feat, describes a creative new application of 
eclipse data, or proposes an ambitious survey to find ever- 
larger numbers of ever-smaller transiting planets. Keeping 
up with the field leaves one breathless, and wary of making 
predictions. 

The only safe bet is that the Kepler space mission will 
have a major impact. In March 2009, the Kepler satellite 
was launched into an an Earth-trailing orbit where it will 
stare at a single field of lO'^ stars for at least 3.5 yr, seeking 
Earth-like planets in Earth-like orbits around Sun-like stars 
(Borucki et al. 2003). Kepler will end our obsession with 
close-in, tidally-influenced, strongly-irradiated giant plan- 
ets. At last it will be possible to study eclipses of multitudes 
of longer-period and smaller-radius planets. Stellar heating 
will not confound models of their atmospheres or interiors, 
and tidal effects will not confound the interpretation of their 
orbital properties. 

Many of Kepler s discoveries will be small enough to 
have solid surfaces, making it possible to measure masses 
and radii of rocky exoplanets. The mission's highest prior- 
ity goal is to determine the abundance of terrestrial planets 
in the "habitable zones" of their parent stars, defined as the 
range of orbital distances within which water could exist 
in liquid form on the surface of a rocky planet. The abun- 
dance of such planets is a key input to the "Drake equation" 
as well as to theories of planet formation and to the designs 
of future missions that will find and study the nearest exam- 
ples. 

The precision of space-based photometry will also be put 
to good use to detect a host of subtle effects that have been 
described in the literature but not yet observed. Already, 
data from the Kepler and CoRoT satellites have been used 
to detect occultations in visible light, as opposed to infrared 
(Borucki et al. 2009, Snellen et al. 2009), although in both 
cases the signal is attributed mainly to thermal emission 
rather than the elusive reflected component. 

Many investigators will seek to detect short-term varia- 
tions in the times, durations, and depths of transits due to 
gravitational perturbations by additional planets, as men- 
tioned in Section 3.2 and discussed further in the chapter by 
Fabrycky. In particular, Kepler might find habitable planets 
not only by detecting their transits, but also by observing 
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Fig. 14. — Examples of data used to measure the projected spin-orbit angle A. The top panels show transit photometry, and the bottom 
panels show the apparent radial velocity of the star, including both orbital motion and the anomalous Doppler shift (the Rossiter- 
McLaughlin effect). The left panels show a well-aligned system and the middle panels show a misaligned system. The right panels 
show a system for which the stellar and orbital "north poles" are nearly antiparallel on the sky, indicating that the planet's orbit is either 
retrograde or polar (depending on the unknown inclination of the stellar rotation axis). References: Winn et al. (2006; 2009a,b). 



their effects on other transiting planets. The same princi- 
ple can be used to detect habitable "exomoons" (Kipping et 
al. 2009) or habitable planets at the Trojan points (Lagrange 
L4 and L5) in the orbits of more massive planets (Ford & 
Gaudi 2006, Madhusudhan & Winn 2009). 

Longer-term orbital perturbations should also be mea- 
surable, due to additional bodies (Miralda-Escude 2002) or 
relativistic precession (Jordan & Bakos 2008, Pal & Koc- 
sis 2008). For very close-in planets, the orbital precession 
rate should be dominated by the effects of the planetary tidal 
bulge, which in turn depends on the deformabiUty of the 
planet. This raises the prospect of using the measured pre- 
cession rate to infer some aspects of the planet's interior 
structure (Ragozzine & Wolf 2009). 

In addition, the precise form of the transit light curve 
depends on the shape of the planet's silhouette. With pre- 
cise enough photometry it may be possible to detect the 
departures from sphericity due to rings (Barnes & Fortney 
2004) or rotation (Seager & Hui 2002, Barnes & Fortney 
2003). Already it has been shown that HD 189733b is less 
oblate than Saturn (Carter & Winn 2009b), as expected if 
the planet's rotation period is synchronized with its 2.2- 
day orbit (slower than Saturn's 1 1 hr rotation period). Soon 
we will have a sample of transiting planets with larger or- 
bital distances, for which synchronization is not expected 
and which may therefore be rotating quickly enough for the 
oblateness to be detectable. 

Another prize that remains to be won is the discovery 
of a system with more than one transiting planet. This 
would give empirical constraints on the mutual inclination 
of exoplanetary orbits (Nesvomy 2009), as well as esti- 



mates of the planetary masses that are independent of the 
stellar mass, through the observable effects of planet-planet 
gravitational interactions. Already there are a few cases in 
which there is evidence for a second planet around a star 
that is known to have a transiting planet (see, e.g., Bakos et 
al. 2009), but in none of those cases is it known whether the 
second planet also transits. 

Ground-based transit surveys wiU stiU play an important 
role by targeting brighter stars over wider fields than CoRoT 
and Kepler. With brighter stars it is easier to rule out false 
positives, and to measure the planetary mass through radial- 
velocity variations of the host star. For the smallest plan- 
ets around the relatively faint stars in the CoRoT and Ke- 
pler fields this will be difficult. Bright stars also offer more 
photons for the high-precision follow-up investigations that 
make eclipses so valuable. 

Ground-based surveys will continue providing such tar- 
gets, and in particular will mine the southern sky, which is 
comparatively unexplored. Two contenders are the Super- 
WASP and HAT-South surveys, which will aim to improve 
their photometric precision and discover many Neptune- 
sized or even smaller planets in addition to the gas giants. 
The MEarth project, mentioned in Section 5.1, is specif- 
ically targeting small and low-luminosity stars because of 
the less stringent requirements on photometric precision 
and because the habitable zones occur at smaller orbital 
distances, making transits more likely and more frequent 
(Nutzman & Charbonneau 2008). 

An attractive idea is to design a survey with the fine 
photometric precision that is possible from space, but that 
would somehow survey the brightest stars on the sky in- 
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stead of those within a narrow field of view. Several mission 
concepts are being studied: the Transiting Exoplanet Sur- 
vey Satellite (Deming et al. 2009) and the PLATO mission 
(Lindberg et al. 2009) would tile the sky with the fields-of- 
view of small cameras, while LEAVITT would have a series 
of telescopes on a spinning platform (OlUng 2007). An- 
other concept is to deploy an armada of "nanosatellites" to 
target individual stars (S. Seager, priv. comm. 2009). 

Each time a terrestrial planet is found transiting a bright 
star, a period of intense anticipation will ensue, as as- 
tronomers try to characterize its atmosphere through trans- 
mission or emission spectroscopy. The anticipation will 
be especially keen for any planets in the habitable zone 
of their parent stars, which may reveal water or other 
molecules considered important for hfe. Attaining the nec- 
essary signal-to-noise ratio will be excruciating work, and 
may require the James Webb Space Telescope or an even 
more capable special-purpose instrument. 

Transits and occupations have had a disproportionate 
impact early in the history of exoplanetary science. This is 
mainly because of the unexpected existence of close-in gi- 
ant planets, and their high transit probabilities. In the com- 
ing years, other techniques such as astrometry, microlens- 
ing, and spatially-resolved imaging will mature and give a 
more complete accounting of exoplanetary systems. Still, 
in closing, let us recall that most of our most fundamental 
and precise information about stars comes from echpsing 
systems, even after more than a century of technological 
development since eclipses were first observed. The same 
is likely to be true for exoplanets. 
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